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Abstract: The IceCube Neutrino Observatory at the South Pole has been completed in December 2010. In this paper we 
describe the final detector and report results on physics and performance using data taken at different stages of the yet 
incomplete detector. No signals for cosmic neutrinos from point sources and diffuse fluxes have been found. Prospects 
of these searches, including the setup of multi-messenger programs, are discussed. The limits on neutrinos from GRBs, 
being far below model predictions, require a reevaluation of GRB model assumptions. Various measurements of cosmic 
ray properties have been obtained from atmospheric muon and neutrino spectra and from air shower measurements; 
these results will have an important impact on model developments. IceCube observed an anisotropy of cosmic rays 
on multiple angular scales, for the first time in the Southern sky. The unique capabilities of IceCube for monitoring 
transient low energy events are briefly discussed. Finally an outlook to planned extensions is given which will improve 
the sensitivities both on the low and high energy side. The IceCube contributions to this conference (ICRC 2011) can be 
found in [1]. 

Keywords: Cosmic neutrinos, cosmic rays, IceCube, DeepCore, IceTop 



1 Introduction 

The main component of the IceCube Neutrino Observatory 
at the geographic South Pole is a 1-km^ detector instru- 
mented with optical sensors in the clear ice of the polar 
glacier at depths between 1450 and 2450 m. The installa- 
tion of IceCube with all its components was completed in 
December 2010. The main purpose of IceCube is the de- 
tection of high energy neutrinos from astrophy sical sources 
via the Cherenkov light of charged particles generated in 
neutrino interactions in the ice or the rock below the ice. 

The basic motivation for the construction of IceCube is to 
contribute to answering the fundamental, still unanswered 
question of the origin of cosmic rays. If cosmic rays are 
accelerated in astronomical objects, like Supernova Rem- 
nants (SNR), Active Galactic Nuclei (AGN) or Ganmia 
Ray Bursts (GRB), one expects the accelerated particles 
to interact with the accelerator environment leading mainly 
to pion production. The principle of such a reaction of an 
accelerated hadron N with an ambient hadron or photon is: 

N+N'n^X+[ ^l^t^^^,^e+^,^,^e (+C.C.) 
[ TT^ ^ 77 

(1) 

While neutral pions decay to gammas which can be de- 
tected by satellite ganmia detectors up to several 100 GeV 
and by Cherenkov ganmia ray telescopes in the TeV range, 
the charged pion decays or other weak decays such as kaon 
decays lead to neutrinos with a similar energy spectrum. 



If the pion production happens in or near the accelera- 
tor one expects to observe neutrino point sources. Inter- 
actions on the interstellar or intergalactic radiation back- 
ground would lead to a diffuse flux of neutrinos. Also the 
summed flux of many faint sources could be seen as diffuse 
flux. The highest energies in the diffuse flux are expected 
to be in the EeV range stemming from interactions of the 
highest energy cosmic rays with the photons of the Cos- 
mic Microwave Background (CMB). The observation of 
these neutrinos could confirm that the cosmic rays are lim- 
ited at energies of about 10^^ eV by the so-called "Greisen- 
Zatsepin-Kuzmin cut-off" (GZK cut-off). The importance 
of the observation of neutrinos from astrophysical sources 
to constrain theoretical models was stressed in various talks 
at this conference, for example [2]. 

In the lowest part of the IceCube detector a subvolume 
called DeepCore is more densely instrumented lowering 
the energy threshold from about 1 TeV in most of the de- 
tector to about 10 GeV. This addition to the original detec- 
tor design extends appreciably the physics reach of the ob- 
servatory to atmospheric neutrino oscillation phenomena, 
WIMP searches at lower masses and improves the sensitiv- 
ity for the detection of transient events like supemovae and 
GRBs. 

IceTop, the surface component of IceCube, is an air shower 
array covering an area of 1 km^ . With this detector air- 
showers from primary particles in the energy range from 
about 300 TeV to above 1 EeV can be measured. The 
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Table 1 : List of the years when a certain configuration of 
IceCube (IC), IceTop (IT) and DeepCore (DC) became op- 
erational. The DC strings are also included in the numbers 
for IC. In this paper we will use abbreviations like IC40, 
IT40 for the constellation in 2008, for example. 



Year 


IC strings 


IT stations 


DC strings 


2006 


9 


9 




2007 


22 


26 




2008 


40 


40 




2009 


59 


59 




2010 


79 


73 


6+7 


2011 


86 


81 


8+12 




IceTop 

'81 Stations, each with 

2 IceTop Cherenkov detector tanks 
2 optical sensors per tank 
324 optical sensors 



IceCube Array 

86 strings including 8 DeepCore strings 
60 optical sensors on each string 
51 60 optical sensors 

December, 2010: Project completed, 86 strings 



DeepCore 

strings-spacing optimized for lower energies 

30 optical sensors 



Figure 1 : The IceCube detector with its components Deep- 
Core and IceTop in the final configuration (January 2011). 



detector is primarily designed to study the mass composi- 
tion of primary cosmic rays in the energy range from about 
10^^ eV to 10^^ eV by exploiting the correlation between 
the shower energy measured in IceTop and the energy de- 
posited by muons in the deep ice, see [3]. 

In the following I will describe the IceCube detector with 
the sub-components DeepCore and IceTop. During the 
construction time from 2004 to the end of 2010 data have 
been taken with the still incomplete detector, see Table 1 . 
Results obtained with differently sized detectors will be re- 
ported for neutrino point source searches, for diffuse neu- 
trino fluxes searches, search for "exotic" (beyond standard 
model) particles and studies of cosmic rays. The summary 
includes a brief outlook to possible extensions in the future. 

2 Detector 



a 10^' photo multiplier tube (PMT) to record the Cherenkov 
light of charged particles traversing the ice. In addition, a 
DOM houses complex electronic circuitry supplying sig- 
nal digitisation, readout, triggering, calibration, data trans- 
fer and various control functions [4] . The most important 
feature of the DOM electronics is the recording of the ana- 
log waveforms in 3.3 ns wide bins for a duration of 422 ns. 
The recording is initiated if a pulse crosses a threshold of 
0.25 photoelectrons. With a coarser binning a 'fast ADC 
extends the time range to 6.4 /is. 

Ice properties: At the depth of the detector the ice is very 
clear with an absorption length reaching up to 200 m. Scat- 
tering and absorption show a depth dependence, which fol- 
lows the dust concentration in the polar glacier. The most 
prominent feature is a significantly higher dust concentra- 
tion at depths between 2000 and 2100 m. The measurement 
and modelling of the ice properties for reconstruction and 
simulation is discussed in [5]. 

DeepCore: In the lower part of the detector a section 
called DeepCore is more densely instrumented. The Deep- 
Core subarray includes 8 (6) densely instrumented strings 
optimized for low energies plus the 12 (7) adjacent standard 
strings (the numbers in parentheses apply to the DeepCore 
configuration of the 2010 running with 79 strings for which 
we will discuss results below). 

IceTop: The 1-km^ IceTop air shower array [3] is located 
above IceCube at a height of 2835 m above sea level, cor- 
responding to an atmospheric depth of about 680 g/cm^. It 
consists of 162 ice Cherenkov tanks, placed at 81 stations 
mostly near the IceCube strings (Fig. 1). In the center of 
the array, three stations have been installed at intermedi- 
ate positions. Together with the neighbouring stations they 
form an in-fill array for denser shower sampling yielding a 
lower energy threshold. Each station comprises two cylin- 
drical tanks, 10 m apart. 

Each tank is equipped with two DOMs to record the 
Cherenkov light of charged particles that penetrate the tank. 
DOMs, electronics and readout scheme are the same as for 
the in-ice detector. The two DOMs in each tank are op- 
erated at different PMT gains to cover linearly a dynamic 
range of about 10^ with a sensitivity to a single photoelec- 
tron). The measured charges are expressed in units of 'ver- 
tical equivalent muons' (VEM) determined by calibrating 
each DOM with muons (see ref. [6]). 



IceCube: The main component of the IceCube Observa- 
tory is an array of 86 strings equipped with 5160 light de- 
tectors in a volume of 1 km^ at a depth between 1450m 
and 2450 m (Fig. 1). The nominal IceCube string spacing 
is 125 m on a hexagonal grid (see DeepCore below). 

Each standard string is equipped with 60 light detectors, 
called 'Digital Optical Modules' (DOMs), each containing 



Trigger and data acquisition: The average noise rate of 
a single DOM is 540 Hz. To initiate the full readout of 
DOMs, a so-called 'hard local coincidence' (HLC) is re- 
quired. In IceCube at least one of the two nearest neigbour 
DOMs of a string must have signals above threshold within 
±1 /iS, resulting in a launch rate of 5-20 Hz per DOM, 
falling with increasing depth. In IceTop the HLC require- 
ment is a coincidence of the two high gain DOMs of a sta- 
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tion. This results in a launch rate of high gain DOMs of 2-4 
Hz compared to about 1600 Hz of a single high gain DOM 
at a threshold of about 0.2 VEM. 

In the counting house at the surface, triggers are formed 
from the HLCs deciding if the data are written to a perma- 
nent storage medium to make it available for later analy- 
sis. The basic in-ice trigger, for example, requires that at 
least 8 DOMs are launched by an HLC leading to a rate of 
about 2 kHz. A very loose trigger requirement is applied to 
the DOMs in the DeepCore fiducial region (below the dust 
layer) by requiring 3 or more HLC hits within a 2.5 /iS time 
window. The basic trigger for IceTop is issued if the read- 
outs of 6 or more DOMs are launched by an HLC leading 
to a rate of 30 to 40 Hz. For all detector components HLC 
hits are always stored in case of a trigger issued by another 
detector component. 

For each DOM above threshold, even without a local co- 
incidence, condensed data, so-called SLC hits ('soft local 
coincidence'), are transmitted. These data contain only the 
time and amplitude of the PMT waveform peaks (for in- 
ice DOMs) or the time and integrated charge (for IceTop 
DOMs). The SLC hits are, for example, used for detect- 
ing transient events and to generate vetos for special event 
signatures. In the case of IceTop they are useful for de- 
tecting single muons in showers where the electromagnetic 
component has been absorbed (low energies, outer region 
of showers, inclined showers). 

For monitoring transient events via rate variations, the time 
of single hits are histogrammed. In IceTop the single hits in 
different tanks are obtained with various thresholds ('scaler 
rates' for heliosperic physics). 

Triggered events which fulfill further criteria for various 
event classes ('muon', 'cascade' etc.) are sent via satellite 
to the IceCube Computing Center in Madison. In addition 
fast online processing produces alerts for other telescopes 
in case of significant neutrino accumulations (see Section 
4.5 on follow-up programs). 



3 Detection Methods and Performance 

IceCube performance: Neutrino point source searches 
rely on directional information and are therefore performed 
with muon neutrinos. The directions of muon tracks can 
be reconstructed to 1° or better (see Moon shadow anal- 
ysis [7]). Primarily IceCube is designed to measure up- 
going neutrinos^ using the Earth as filter against the large 
background of high energy muons from cosmic rays. How- 
ever, because the neutrino cross section increases with en- 
ergy the Earth becomes opaque for neutrinos above about 
1 PeV. This can be seen in Fig. 2 where the energy depen- 
dence of the effective area of IceCube is plotted for differ- 
ent zenith angles. The effective area is defined as the target 
area which yields the observed muon neutrino rate when 
each neutrino is detected with 100% probability. Since at 
high energies the background from down-going cosmic ray 
muons becomes relatively small IceCube has extended its 
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Figure 2: Effective area for muon neutrino detection with 
IceCube (IC40) as a function of the neutrino energy for dif- 
ferent zenith angle ranges of the Northern sky. The plot was 
made for the IC40 analysis of diffuse neutrino fluxes [8]. 



search also to the Southern sky at high energies (see Sec- 
tion 4). 

For up-going neutrinos the background is mainly atmo- 
spheric neutrinos generated in the Northern atmosphere, 
while the background for down-going neutrinos comes 
mostly from high energy atmospheric muons. The extrac- 
tion of signals for astrophysical neutrinos relies either on 
accumulations in space (point sources or extended objects) 
and/or in time (flares, GRB etc.) or on the assumption that 
the cosmic neutrino spectra are harder than for secondary 
cosmic rays, often a spectral index of about -2 compared to 
-3.7 for atmospheric muons and neutrinos is assumed. The 
latter is particularly important for the measurement of dif- 
fuse neutrino fluxes where an on/off source determination 
of background is not possible. 

The muon neutrino energy cannot be directly measured 
(except in the cases where a neutrino interacts in the de- 
tector and the muon ranges out). The measured energy 
loss of muons in the ice is used as a proxy for the neu- 
trino energy. For muon energies above about 1 TeV, cor- 
responding to the critical energy of muons in ice, the to- 
tal energy loss rises approximately linearly with the muon 
energy due to bremsstrahlung, pair production and nuclear 
interactions. This allows the determination of the muon en- 
ergy from the energy loss with a resolution of about 50% 
(Alog^Q E ^ 0.2). The muon energy yields only a very 
coarse proxy for the neutrino energy which is only partially 
transferred to the muon. The angular resolution for muon 
neutrinos is about 1° at 1 TeV and about 0.5° at 1 PeV. 

All flavours of neutrinos can contribute to the search of dif- 
fuse fluxes if they generate an electromagnetic or hadronic 
cascade in the ice. Electron and tau neutrinos can generate 
electromagnetic cascades in charged current interactions, 
all flavours can generate hadronic cascades via neutral cur- 
rent and charged current interactions. Cascades appear in 
IceCube as nearly spherical isotropic light sources, so that 
little direction information can be obtained. On the other 
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Figure 3: Skymap of neutrino candidates in equatorial co- 
ordinates for the IC40+59 data sample (left). The curved 
line indicates the galactic plane. The declination depen- 
dence of the selected energy ranges is shown on the right. 



hand, the neutrino energy resolution is much better than 
in the case of neutrino detection via muons, about 30% at 
10 PeV (Alogio^^ ^ 0-13). Neutrinos have to interact 
in or near the detector to be detectable as cascades which 
makes the neutrino effective area about an order of magni- 
tude lower than for muon tracks. 

Deep Core performance: The main improvement added 
by DeepCore is the decrease of the energy threshold to 
about 10 GeV. Using the surrounding IC strings as veto, 
one can identify low energy neutrino interactions inside or 
near the DC volume. In this way, atmospheric neutrinos 
can be collected over an angular range of 47r sr, yielding un- 
precedented statistical samples of atmospheric neutrinos, 
about 150,000 triggered atmospheric muon neutrinos per 
year, thus allowing oscillation studies [9] . The observation 
of a sizeable number of cascade event in DeepCore [10] 
confirms the performance expectations. 

IceTop performance: IceTop will cover a primary en- 
ergy range from about 300 TeV to 3 EeV for zenith angles 
up to about 65° . In coincidence with IceCube the zenith an- 
gle range is more limited yielding an angular coverage of 
about 0.3 sr. The event rate is sufficient for a composition 
analysis up to about 1 EeV. 

The following resolutions have been obtained for 10 PeV 
(100 PeV) and for zenith angles smaller than 30° [11]: core 
position 7 m (8 m), zenith angle 0.5° (0.3°), energy resolu- 
tion 0.05 (0.04) for log^o E/VeV. 

4 Neutrino Point Sources 
4.1 Search Strategy 

The neutrino point source search relies on the precise di- 
rectional information from muons generated by muon neu- 
trinos interacting in the ice in and around the detector or the 
Earth crust below the detector. Figure 3 shows a skymap of 
arrival directions of neutrino candidates. The plot contains 
57460 up-going and 87009 down-going neutrino candi- 
dates selected from 723 days of data taken with the 40 and 



59 string configurations during 2008 and 2009 (IC40+59). 
In the Northern sky high energies are limited by neutrino 
absorption in the Earth (Fig. 2), in the Southern sky the en- 
ergy threshold has to be increased to reject the large back- 
ground from atmospheric muons. 

In an unbiased search each direction has to be scanned lead- 
ing to a large number of trials and thus a significance reduc- 
tion. To improve signal significances one wants to reduce 
the number of trials by using additional information on the 
signal probabilities: 

- Predefine a list of candidate sources which are theoreti- 
cally likely to emit neutrinos. 

- The list search can be further improved by summing the 
fitted signals for many sources ('stacking'). 

- Search for extended sources on scales from a few de- 
grees, just resolvable, to scales of the size of the galactic 
plane. 

- Search for time and spatial correlation with transient 
events, like flares in AGN. 

- A special class of transient events are GRBs with par- 
ticularly short emission times. The similar properties of 
different GRBs make them well suited for stacking (Sec- 
tion 4.4). 

- Since IceCube is sensitive about 99% of the time to the 
full sky, alerts can be given to other telescopes if IceCube 
detects multiplets of neutrino candidates which accumu- 
late in space and time. Such 'Follow-up Programs' are 
realised with optical. X-ray and 7-ray telescopes. 

4.2 Full Sky Time-Integrated Search 

In a basic approach one searches in the full considered data 
set for a significant accumulation of events in an angular 
range compatible with the angular resolution. For that pur- 
pose a likelihood function is defined which takes into ac- 
count a possible signal and background: 

N 

iK,7)=n[|^^+(i-|)^^] 

For a given direction on the sky Si and Bi are the proba- 
bilities for event i to be signal or background, respectively; 
N is the number of events which is looped over and Us 
is the number of most likely signal events. The likelihood 
function depends also on the energy via a spectral index 7 
which is estimated in the search procedure. The search has 
to be done in a fine grid of directions, here at about 100 000 
points, which reduces appreciably the "pre-trial" signifi- 
cance for a point source to a "post-trial" significance. The 
significances are evaluated defining a test statistics which 
compares the most likely values fis , 7 with the null hypoth- 
esis. Using simulations the distribution of the test statistics 
for the case of no signal is evaluated yielding a p-value 
which is the probability to reach the observed or a higher 
significance for a result , 7 if there is no signal. 

In the analysis of the IC40+59 data the hottest spot at 
(Ra, Dec) = (75.45°, -18.15°) has a pre-trial p-value of 
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Figure 4: Neutrino point source limits (90% c.l.) for an 
E~'^ spectrum. The currently most stringent limit from 
IC40-F59 data is compared to previous and expected lim- 
its. 

Ppre = 10~^'^^, corresponding to an about 4 sigma signif- 
icance, but a post- trial p- value of 0.67, indicating a high 
compatibility with the null hypothesis. This means that no 
significant point source observation can be reported from 
this search. 

An overview of limits obtained from time integrated point 
source searches is given in Fig. 4. The IceCube 40-1-59 re- 
sults are compared to previously published limits from Ice- 
Cube and other experiments. The recently published IC40 
results [12] include also limits for specific source candi- 
dates which had been selected before looking at the data 
(see the list and more details in [12]). With these IceCube 
measurements the limits decreased by about a factor 1000 
over the last 15 years. 

The IC40-F59 limits reached already the projected sensi- 
tivities obtainable by the full detector in one year. How- 
ever, sensitivities below about 10~^^£^~^TeV~^cm~^s~^ 
are necessary to seriously scrutinize models for cosmic ray 
acceleration with neutrino production. Hence in this search 
mode several years of additional data taking might be nec- 
essary for either a neutrino source discovery or a rejection 
of the models. 

4.3 Time-Dependent Searches for Point Sources 

The statistical significance can be improved by includ- 
ing time dependence in the likelihood function (2) since 
sources such as AGN can exhibit significant time variabil- 
ity in photon fluxes, which might also be visible in neutri- 
nos, while the atmospheric background is roughly constant. 

An example for an 'untriggered' search, i.e. without a pri- 
ori time information, was presented at this conference [13] 
using the IC40-F59 data as for the time-independent search 
(the IC40 analysis has recently been published [14]). The 
time-dependent likelihood term for this search is a Gaus- 
sian function, with its mean and width as free parameters. 
Scanning for flares of all durations from 20 seconds to 150 
days the likelihood maximization returns the most signifi- 
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Figure 5: IC59: The time distribution of the signal-to- 
background ratio of events from the location of maximum 
significance. The curve is the fitted gaussian for the most 
significant flare. 



cant flare from a particular direction. The strongest devia- 
tion from background was found in the IC59 data in a di- 
rection (Ra, Dec) = (21.35°, -0.25°), centered on March 
4, 2010 with a FWHM of 13 days (Fig. 5). An excess of 
14 events is seen with a soft spectrum of E~^-^, i.e. with 
no discrimination against the atmospheric spectrum. The 
post-trial p-value is determined to be 1.4% (corresponding 
to about 2.3 sigma) which is not sufficiently significant for 
claiming a neutrino flare discovery. 

4.4 Gamma Ray Bursts 

The large energy dissipation in a Gamma Ray Burst (GRB) 
of about 10^"^ J suggests that a large fraction of the extra- 
galactic cosmic rays at the highest energies could be ac- 
celerated in GRBs. GRBs are usually modeled as explo- 
sions of very massive stars which eventually collapse to a 
black hole. In such models the observed gamma rays stem 
from synchrotron radiation and/or inverse Compton scat- 
tering of electrons accelerated in shock fronts in the colli- 
mated explosive outflow. It was proposed that in the same 
way also protons are accelerated [15, 16]. These protons 
would undergo interactions with the surrounding photon 
field in the fireball and thus generate neutrinos according 
to (1). With their preferred parameters the models predict 
that GRB neutrinos be detectable by IceCube in less than 
one year. 

At this conference a search using IC59 data was reported 
[17]. The search was based on a list of 98 GRB observa- 
tions reported from satellites during times when IceCube 
was taking data. The neutrino search was done similarly to 
the point source searches, using a likelihood like (2) with an 
additional term for the time included. The time probability 
density function was flat in the interval were the first and 
last gamma rays were observed falling off smoothly to both 
sides. In the point- spread function the uncertainty in the 
GRB coordinates as obtained from satellites was included. 
No neutrino candidate was observed in the space-time win- 
dows. The analysis sets a limit far below the predicted 
model flux (Fig. 6). Combining the results from IC40 [18] 
and IC59 the corresponding upper limit lies a factor 5 be- 
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Figure 6: Limits on neutrino flux from GRBs compared to 
models from Waxmann [15] and Guetta et al. [16]. The 
derivation of the limits is based on the Guetta et al. model 
and accounts for the estimated properties of individual 
GRBs (the Waxmann predictions use average properties). 
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Figure 7: Limits and predictions for diffuse muon neutrino 
fluxes. 



low the model curve. This leads to the conclusion that 
either the model picture of GRBs is wrong or the chosen 
parameter values are not correct. Important model param- 
eters are the Lorentz boost factor F of the collimated out- 
flow of the exploding star and the typical time scale t^^^ of 
subsequent collisons of internal shocks. In [17] the limits 
obtained for the combination of these parameters are pre- 
sented. 

4.5 Follow-Up Programs 

A special feature of the IceCube detector is that it is able to 
monitor the whole sky (though with different energy sen- 
sitivities, see Section 3). This can be exploited to send 
alerts to other telescopes with narrow fields of view (op- 
tical, X-ray, gamma-ray) if in a certain space-time window 
an excess of neutrinos above background is observed with 
a predefined significance. The alerted telescopes can then 
make follow-up observations on these 'targets of opportu- 
nity' which would lead to a significance enhancement if 
a positive correlation between different messenger signals 
are observed. The alert decisions have to be made fast, 
i.e. online at Pole and reported via satellite, and have to be 
tuned in a way that the alert rate is tolerable for the alerted 
partners. 

The IceCube collaboration has follow-up programs estab- 
lished with several telescopes: 

- Search for GRB and core-collapse supernovae: neu- 
trino multiplets in a short time window, < 100 s, gen- 
erate alerts for optical follow-up by the Robotic Optical 
Transient Search Experiment (ROTSE) and the Palomar 
Transient Factory (PTF), see [19]. Furthermore an X- 
ray follow-up by the Swift satellite of the most signifi- 
cant multiplets has been set up and started operations in 
February 2011 [20]. 

- Search for neutrinos from TeV-gamma flares: a follow- 
up program with the MAGIC telescope has been tested 
with the IC79 setup and should become active for the 
IC86 running [21]. 



5 Diffuse Flux of Neutrinos 

If there are many point sources, each with an unobservably 
low flux, then the aggregate flux may still be observable 
as a diffuse flux. Interactions of the cosmic rays with the 
matter and radiation near the source or somewhere else on 
their path through the space would lead, according to eq. 
(1), to meson production and the subsequent weak decays 
to a diffuse flux of neutrinos. 

The identification of diffuse cosmic neutrinos relies on the 
assumption that they have a harder spectrum, e.g. E~'^ 
compared to about E~^-^ for atmospheric neutrinos. The 
'prompt' component of atmospheric neutrinos from decays 
of heavy flavour hadrons, which are produced predomi- 
nantly in the first interactions in the atmosphere, is pre- 
dicted to be harder than the 'conventional' neutrino flux. 
This introduces some additional uncertainty in the transi- 
tion region where the cosmic flux is expected to become 
dominant. The experimental limits tell us that this transi- 
tion is well above 100 TeV neutrino energy (see Fig. 7). 

5.1 Diffuse Muon Neutrino Flux 

Figure 7 shows the currently best limit obtained from IC40 
data of up-going muons [22]. The points are the atmo- 
spheric neutrino spectrum determined by unfolding the 
measured muon energy depositions to obtain the flux as a 
function of the neutrino energy. The limit is now below 
the Waxmann-Bahcall bound [24] which gives a guideline 
of how much flux can be at most expected if cosmic neutri- 
nos are generated in or near the accelerating sources (AGN, 
GRB etc.) via meson production as in (1). 

5.2 Cascades and All-Flavour Neutrino Flux 

The interaction of electron neutrinos in IceCube generates 
an electromagnetic cascade which manifests itself in the 
detector as a nearly spherical source of light with little 
information about the direction. To this 'cascade chan- 
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Figure 8: All-flavor flux limit from IC40 data derived from 
a cascade analysis [25] (blue line) and from an EHE neu- 
trino search [27, 28] (red lines). The model dependent lim- 
its assuming an neutrino spectrum are shown as hori- 
zontal lines. The curved red line is the model independent 
quasi-differential limit (normalized by decade in energy). 
The EHE neutrino flux limit is shown together with lim- 
its of other experiments (employing radio techniques). The 
estimated reach for the full detector in 1 and 5 years (black 
dashed lines) is compared to the specific model [23] shown 
as band on the plot. 



cosmic rays are limited at energies of about 10^^ eV by 
the "GZK cut-off", at the point where the 7cmb - nu- 
cleon system surpasses the threshold for pion production 
(with a strong enhancement due to the A -resonance close 
to threshold). Since all involved processes and particles 
are well known this GZK process could be considered a 
'guaranteed' source of cosmogenic neutrinos. However, 
in detail the theoretical predictions for the fluxes vary by 
about 3 orders of magnitude, depending mostly on the as- 
sumed primary composition and the distribution of cosmic 
ray sources. 

At this conference results for 'Extremely-High Energy' 
(EHE) neutrinos using the IC40 detector have been pre- 
sented [27, 28]. The analysis aims at finding down-going 
neutrinos generating very bright events in the detector. 
However, the large atmospheric muon background restricts 
the search to events coming from near the horizon where 
neutrinos have also the largest interaction probability. The 
obtained EHE neutrino flux limits are shown in Fig. 8. Up 
to about 10 EeV IceCube has the best upper limits. The 
comparison with predictions shows that a positive obser- 
vation of GZK neutrinos might still take some years. On 
the other hand improvements in the analysis procedure 
could increase the detection efficiency [29]. For example 
a scheme is currently investigated to use single-tank hits 
in IceTop for a veto against the overwhelming background 
from downgoing muons [30] . 



nel' also neutral current interactions of all flavours, gen- 
erating hadronic cascades, contribute. Therefore the re- 
sults for cosmic neutrinos are expressed as all-flavour neu- 
trino fluxes assuming a flavour ratio of 1 : 1 : 1 at the detector 
(evolving from a 1:2:0 ratio at the source by mixing). For 
diffuse flux measurements, the inferior directional resolu- 
tion is not a major draw-back, but the relatively good en- 
ergy resolution has substantial advantages (see Section 3). 

At this conference a result on a cascade analysis using IC40 
data taken over 374 days was presented [25]. Above a cas- 
cade energy cut of 25 TeV, 14 events are left in data while 
atmospheric neutrino simulations predict less than 4 events. 
Visual inspections of the experimental events indicate that 
at least a fraction of the excess events are induced by at- 
mospheric muons. Because of insufficient statistics in the 
atmospheric muon background simulations and because of 
possible inaccuracies of the atmospheric neutrino flux pre- 
diction, this result is, at this moment, inconclusive (see also 
[26]). 

The current best limit from cascades as derived in a more 
conservative analysis of the IC40 data for an energy range 
between about 90 TeV and 20 PeV is shown in Fig. 8. 

5.3 Extremely-High Energy Neutrinos 

Interactions of the highest energy cosmic rays with the pho- 
tons of the Cosmic Microwave Background (CMB) are pre- 
dicted to generate a diffuse neutrino flux in the EeV range. 
The observation of these neutrinos could confirm that the 



6 Exotics 

An essential part of the IceCube physics program deals 
with generic Particle Physics problems such as the search 
for new particles beyond the Standard Model, called 'ex- 
otic particles'. Such particles include Dark Matter candi- 
dates such as proposed by Supersymmetry (SUSY) or by 
Kaluza-Klein models. The breaking of larger symmetries, 
as postulated by 'Grand Unified Theories' , implies the gen- 
eration of topological defects such as monopoles which can 
also be searched for with IceCube. 

6.1 WIMP Search 

It is now experimentally well established that 'Dark Mat- 
ter' (DM) exceeds normal, baryonic matter by about a 
factor 5. In most common scenarios the DM consists of 
Weakly Interacting Massive Particles (WIMPs) which re- 
mained from the Big Bang after the expansion rate of the 
Universe surpassed their annihilation rate. A promising 
WIMP candidate is the lightest supersymmetric particle, in 
most SUSY variants the neutralino x- In the searches re- 
ported below parameters within the MSSM ('Minimal Su- 
persymmetric Model') have been investigated. 

There are three general DM search strategies: In direct 
searches one looks for elastic WIMP scattering off nuclei; 
in indirect searches one tries to detect WIMP annnihila- 
tion products, such as gammas or neutrinos, with astropar- 
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Figure 9: Preliminary limits on the WIMP induced muon 
flux from the Sun. See text for more explanations. 

tide detectors and finally in accelerator experiments one 
searches for pair-production of DM candidates. None of 
these searches has yet resulted in a conclusive discovery. 

WIMPs from the Sun: IceCube is looking for neutri- 
nos from WIMP (or other exotics) annihilation. The as- 
sumption is that WIMPs would accumulate in gravitational 
wells like the Sun or the center of the Milky Way. In 
Fig. 9 limits on an excess flux of muons from the Sun are 
given for WIMP masses from 50 GeV to 5 TeV [31]. The 
excess determination assumes a muon neutrino spectrum 
which depends on the WIMP mass and the annihilation 
channel. The studied channels VK+l^" and bb have par- 
ticularly hard and soft spectra, respectively (the harder the 
spectrum, the higher IceCube 's sensitivity). The analysis 
combines data taken with IceCube and the precursor detec- 
tor AMANDA between 2001 and 2008, with a total live- 
time of 1065 days when the Sun was below the horizon 
[31]. In 2008 AMANDA was taking data together with 
IceCube in the IC40 configuration employing an integrated 
data acquisition system. The figure also shows the esti- 
mated sensitivity for the full detector. 

The muon flux limits can be related to direct measure- 
ments using the following chain of arguments: The ac- 
cumulation requires a cool-down of the WIMPs by elas- 
tic scatters to get gravitationally trapped. The capture rate 
for WIMPs with a fixed mass can be calculated for given 
elastic cross sections (for different nuclei, with and with- 
out spin dependence) and the local WIMP density and ve- 
locity distribution. In this analysis the 'canonical' values 
Pbm = 0.3GeV/cm^ and ^dm = 270km/s have been 
asummed. Assuming further that an equilibrium between 
capture and annihilation is reached in the Sun, the annihi- 
lation rate becomes half of the capture rate per WIMP. With 
further assumptions about the annihilation channels the rate 
of muons from muon neutrinos can be determined as de- 
scribed above. The shaded area in Fig. 9 indicates the re- 
gion not yet excluded by the MSSM parameter constraints 
through the direct searches by the experiments CDMS and 
XENONIOO. For more details see [31]. 



WIMP annihilation in the Milky Way and dwarf 
spheroidals: In the contribution [32] a search for a neu- 
trino excess from the galactic center and halo was reported. 
Using IC40 data (367 days) preliminary limits on {aann 
have been obtained in the range 10~^^ to 10~^^ cm^ s~^, 
depending on the WIMP mass. The 'natural scale', derived 
from cosmological parameters, is about 3 • 10~^^ cm^ s~^. 
The limits depend strongly on the assumed model for the 
WIMP density and on the annihilation channel. A compar- 
ison of the limits for the r+r~ channel with the regions 
preferred by the satellite experiments PAMELA and Fermi 
(see details in [32]) shows that the WIMP searches of Ice- 
Cube are constraining these preferred regions. 

Another WIMP search reported at this conference is aiming 
at spheroidal dwarf galaxies [33] yielding the sensitivities 
for IC59 data. Although the sensitivities do not yet reach 
those of the 7-ray measurements (MAGIC, Fermi) the neu- 
trino channel certainly adds complementary information. 

6.2 Magnetic Monopoles 

Relativistic magnetic monopoles, if they exceed the 
Cherenkov threshold at /3 0.76, deposit huge amounts of 
light in the detector and thus have a very clear signature. In 
a contribution to the conference a preliminary upper limit 
for the monopole flux was reported using IC22 data [34] . 
This limit is orders of magnitude better than previous ones. 

Also discussed in [34] are the prospects of improving this 
limit with IC40 data to a level which is about a factor 1000 
below the Parker bound. The Parker bound relates the ob- 
served strengths of cosmic magnetic fields to the maximal 
possible abundance of monopoles exploiting that the accel- 
eration of monopoles by magnetic fields would damp those 
fields. The non-observation of monopoles constrains the 
combinations of Grand Unified Theories and inflationary 
scenarios. 

7 Cosmic Rays 

Origin, composition and spectrum of high energy cosmic 
rays are still not well understood. In particular above some 
100 TeV, up to which direct measurements with balloons 
and satellitites are posssible, the experimental situation is 
far from being satisfactory. The understanding of the muon 
and neutrino fluxes from cosmic ray initiated air showers is 
also essential for IceCube because they are the major back- 
ground in the search for extraterrestial neutrinos and exotic 
particles. 

The IceCube observatory offers a variety of possibilities 
to measure cosmic rays, analyse composition and deter- 
mine the spectra which can be used to tune the mod- 
els. IceCube can be regarded as a cubic-kilometer scale 
three-dimensional cosmic ray detector with the air show- 
ers (mainly the electromagnetic component) measured by 
the surface detector IceTop and the high-energetic muons 
and neutrinos measured in the ice. In particular the mea- 
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Figure 10: Demonstration of the composition sensitivity of 
the in-ice muon spectrum by measurements of the muon 
spectrum using the full year of ICS 9 data. Plotted is the ra- 
tio of data over the simulation prediction (details in [35]) 
as a function of the reconstructed muon surface energy, 
which is about a factor 10 lower than the primary cosmic 
ray energy. The comparison indicates a preference for ei- 
ther no change in the spectral slope around the knee ("no 
knee", not favoured by other observations) or the same 
slope change for all elements. 



surement of the dominantly electromagnetic component of 
the airshower in IceTop in coincidence with the high en- 
ergy muon bundle (muon threshold about 500 GeV), orig- 
inating from the first interactions in the atmosphere, has 
a strong sensitivity to composition. Here IceCube offers 
the unique possibility to clarify the cosmic ray composi- 
tion and spectrum in the range between about 300 TeV and 
1 EeV, including the 'knee' region and a possible transition 
from galactic to extra-galactic origin of cosmic rays. 

7.1 Cosmic Ray Physics with Muons in IceCube 

Atmospheric muon spectra: Atmospheric muon and 
neutrino spectra measured with IceCube probe shower de- 
velopment of cosmic rays with primary energies above 
about 10 TeV. In a contribution to the conference [35] it 
was shown that with an accurate measurement of the muon 
spectra one can discriminate between different composition 
models (Fig. 10). At the current stage of the investigation 
a smoother transition of the different element contributions 
in the knee region (than for example suggested by the poly- 
gonato model [36]) is preferred. With additional systematic 
studies a clarification should be reached about what energy 
dependence of composition has to be used in simulation 
models. 

This is a completely new approach to analyse cosmic ray 
composition in the knee region which is otherwise difficult 
to tackle. For the analysis new methods had to be devel- 
oped, for example a method for determination of the energy 
of the leading muon by exploiting cascade-like stochastic 
energy losses [35]. 



Laterally separated muons: At high energies the muons 
reach the in-ice detector in bundles which are, for primaries 
above about 1 PeV, collimated within radii of the order of 
some 10 m. Most of the muons stem from the soft pe- 
ripheral collisons with little transverse momentum transfer. 
Perturbative QCD calculations, however, predict the occur- 
rance of muons with higher transverse momenta in some 
fraction of the events. A first analysis of the IC22/IT26 
data [37], where the muon bundle was measured together 
with the shower energy in IceTop, demonstrated that sepa- 
rations of single muons from the bundle by more than about 
100 m, corresponding to transverse momenta above about 
6 GeV, could be detected. A better understanding of the re- 
maining background from uncorrelated multiple events and 
an unfolding from the lateral separation to transverse mo- 
mentum distributions is currently pursued. With a larger 
detector and also without requiring IceTop coverage, the 
statistics will be sufficient do perform a detailed analysis 
and comparison to the model predictions for meson produc- 
tion. This will have important implications for air shower 
simulations which the cosmic ray analyses have to rely on. 

Seasonal variations of the muon rate: IceCube ob- 
serves a ±8% seasonal variation of muon rates in the ice. 
This modulation is strongly correlated with the variability 
of the temperature, and thus of the density, in the upper 
atmosphere at heights corresponding to pressures around 
10 to 100 hPa. The convolution of the density profile with 
the production cross section for muons defines the effec- 
tive temperature Teff. The relation between the effective 
temperature change and the rate change, assumed to be lin- 
ear. 



aT- 



ATeff 



(3) 



depends on the K/tt production ratio. From the coefficient 
aT measured over 4 years on a sample of 150 billion events 
a preliminary result is reported in [38] which indicates that 
the currently assumed K/tt = 0.15 has to be lowered to 
about 0.1. If confirmed, this would lead to modifications of 
the models for air shower simulation. 

7.2 Atmospheric Neutrino Spectra 

Muon neutrinos: IceCube has the most precise deter- 
mination of the atmospheric muon neutrino spectrum at 
high energies (Fig. 7). This spectrum has to be unfolded 
from the measured muon energies to the neutrino energies. 
A measurement of electromagnetic and hadronic cascades 
generated by neutrinos of all flavours would yield a better 
energy determination. This is important especially at the 
high energy end where signals from diffuse neutrino fluxes 
are searched for. The potential transition from atmospheric 
to cosmic neutrinos is still theoretically uncertain due to 
missing information about composition and about the un- 
certainty in the prompt contribution from heavy quark pro- 
duction. 
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Figure 1 1 : Relative intensity map in equatorial coordinates 
for cosmic rays of the 20-TeV sample (top) and the 400- 
TeV sample (middle). The curve indicates the galactic 
plane with the dot as the galactic center. The projections 
unto the right ascension of both maps are shown at the bot- 
tom. 

Cascade analysis with DeepCore: In a first analysis of 
data taken with the DeepCore detector in the IC79 configu- 
ration (2010/1 1) cascades from atmospheric neutrinos have 
been analysed. In the DeepCore detector 1029 cascade can- 
didates have been observed with a medium energy around 
180 GeV for 281 days of data [10] while 1104 were pre- 
dicted from simulations using the Bartol model [39]. Of 
the predicted events 59% are cascades with about equal 
amounts of Ue CC and NC events. The remaining 41% 
is background from muon tracks from up-going u^. Back- 
ground from down-going atmospheric neutrinos and other 
systematic uncertainties are still under investigation [10]. 

This result from the newly commissionned DeepCore de- 
tector supports the expectations for the performance of the 
detector (Section 3). The physics goals of measuring neu- 
trino oscillations [9], decreasing the mass range for the 
WIMP search and enhancing the sensitivity for supernovae 
detection become very realistic. 

7.3 Cosmic Ray Anisotropy 

IceCube has collected a large amount of cosmic ray muon 
events, about 10^^ events between 2007 and 2010, and ev- 
ery year of running with the full detector will increase this 
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Figure 12: Cosmic ray anisotropics on the scale of 10 to 
30°. Representation as for the maps in Fig. 11 

number by about the same amount. These events have been 
used to study cosmic ray anisotropics, for the first time in 
the Southern sky. The observation of anisotropics on mul- 
tiple angular scales has been previously reported [40, 41]. 
At this conference, analyses of anisotropics using 33-10^ 
events from IC59 data were presented with results on en- 
ergy and angular scale dependencies as well as various sta- 
bility tests of the analyses [42, 43, 44]. 

Figure 1 1 shows skymaps of relative intensities for selec- 
tions of muon energies resulting in primary energy distri- 
butions which center around 20 TeV and 400 TeV. In the 
20-TeV right ascension projection a clear structure dom- 
inated by a dipole and quadrupole contribution is visible 
while the most significant feature in the 400-TeV data set 
is a deep deficite with a completely different phase than the 
dip in the 20-TeV data. For more details see [42] . 

In addition to large-scale features in the form of strong 
dipole and quadrupole moments, the data include several 
localized regions of excess and deficit on scales between 
10° and 30° (Fig. 12). Angular decomposition into speri- 
cal harmonics exhibits significant contributions up to 1=15. 
More details can be found in [43]. 

As yet the anisotropics observed on multiple angular scales 
and at different energies have not found an explanation. 
One could expect an effect due to the movement of the solar 
system relative to the Milky Way, the so-called Compton- 
Getting effect. This effect which results in a dipole compo- 
nent in the cosmic ray intensity distribution cannot, at least 
not fully, explain the data. Theoretical explanations like lo- 
cal magnet fields affecting the cosmic ray streams and/or 
nearby sources of cosmic rays are discussed. The deter- 
mination of the energy dependence of anisotropics will be 
crucial for scrutinizing models. For this reason an analysis 
using IceTop with a better energy resolution and an exten- 
sion to the PeV range for the primary cosmic rays has been 
started. 

7.4 Cosmic Ray Composition 

As mentioned above, the combination of the in-ice detector 
with the surface detector offers a unique possibility to de- 
termine the spectrum and mass composition of cosmic rays 
from about 300 TeV to 1 EeV 

The first analysis exploiting the IceTop-IceCube correla- 
tion was done on a small data set corresponding to only 
one month of data taken with about a quarter of the final 
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Figure 13: Simulated correlation between the energy loss 
of the muon bundels in the ice (K70) and the shower size at 
the surface (S125) for proton and iron showers. The shad- 
ing indicates the percentage of protons over the sum of pro- 
tons and iron in a bin. The lines of constant primary energy 
are labelled with the logarithms of the energies. 
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Figure 14: Average logarithmic mass of primary cosmic 
rays measured with IC40/IT40. 

detector [45]. The energy was restricted to 1 to 30 PeV. A 
neural network was employed to determine from the mea- 
sured input variables shower size and muon energy loss the 
primary energy and mass (Fig. 13). The resulting average 
logarithmic mass is shown in Fig. 14. These results are still 
dominated by systematic uncertainties, such as the energy 
scale of the muons in IceCube and of the effects of snow 
accumulation on the IceTop tanks. 

A first look into the IC79/IT73 data set taken in 2010 shows 
that there will be enough statistics for composition analy- 
sis up to at least 1 EeV [46]. An estimation yields about 
150 events with energies larger than 300 PeV and 15 events 
larger than 1 EeV in 1 year of data taking with the full de- 
tector. 

The systematic uncertainties related to the models can be 
reduced by including different mass sensitive variables, 
like zenith angle dependence of shower size [11], muon 
rates in the surface detector and shower shape variables, 
and checking for consistency. 



The IceCube-IceTop combination has also been used to 
identify high-energy photons as IceTop showers with no 
muons in the ice [47] 

8 Transient Rate IMonitoring 

Transient events such as supernovae, GRBs or sun flares, 
if they generate very high fluxes of low energy particles, 
could be observed as general rate increases above the noise 
level in the DOMs even if they could not be detected indi- 
vidually by IceCube or IceTop. 

Supernova explosions in our and nearby galaxies would be 
observable by IceCube via a rate increase in all DOMs due 
to a high interaction rate of low energy neutrinos. With a 
rather low average noise of 286 Hz per DOM (with after- 
pulse rejection) IceCube is particularly suited to emit su- 
pernova alerts, specifically important when the supernova 
is obscured by dust or stars in a dense region. Measure- 
ments would be sensitive to the supernova parameters such 
as the progenitor star mass, neutrino oscillations and hier- 
archy. In the contribution [48] possibilities for improving 
the current sensitivities, including also DeepCore, are dis- 
cussed. 

IceTop is able to monitor cosmic ray products from tran- 
sient events such as from Sun flares, as demonstrated with 
the observation of the Dec 13, 2006 Sun flare event [49]. 
The detector readout has since then been setup such that 
counting rates could be obtained at different thresholds al- 
lowing to unfold cosmic ray spectra during a flare. At this 
conference the observation of a Forbush decrease in Febru- 
ary 2011 was reported [50]. 

9 Summary and Outlook 

The IceCube Neutrino Observatory has been completed 
and reached or exceeded its design sensitivity. As yet, re- 
sults from the partly completed detector (IC22,40,59) show 
no evidence for cosmic neutrinos, although the detector 
reached sensitivities which are either close to model predic- 
tions or are sometimes seriously challenging models. Most 
notable is the IC40+59 limit on GRBs which is 5 times 
below the model prediction of [16] with preferred parame- 
ters, demanding a reassessment of the model and/or param- 
eters. Point source searches, time dependent or not, with 
and without candidate lists, have not yet reached the level to 
constrain the most common models, but will in some years 
of running. The hope is to accelerate the progress by fur- 
ther developing methods to enhance significances, for ex- 
ample by employing multi-messenger methods and follow- 
up programs with optical. X-ray and 7-ray telescopes. 

The limits on diffuse cosmic neutrino fluxes are now a fac- 
tor of 2.5 below the Waxman-Bahcall bound, indicating 
that the limits have reached a relevant region of predic- 
tions. The reported analyses of cascade events promise 
the opening of a new window for studies of atmospheric 
neutrinos, in particular their 'prompt' contributions, and of 
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cosmic neutrinos with good energy resolution. In the EHE 
region the sensitivity to the range of GZK predictions will 
be reached within a few years. 

Concerning searches for 'exotic' particles, limits for WIMP 
masses between 50 GeV and 5 TeV have reached regions 
in the parameter space which are not excluded by direct 
search experiments. Magnetic monopole limits are now 
nearly a factor 1000 below the 'Parker Bound' (upper 
bound derived from the strength of existing cosmic mag- 
netic fields) and are constraining GUT models. 

Although most of these limits are very important and 
unique complements to results with other messengers it 
is comforting to know that also positive observations have 
been made with IceCube. These results concern mainly the 
field of cosmic rays and are mostly of high importance for 
the improvement of cosmic ray and airshower models. Re- 
sults have been reported on atmospheric neutrino and muon 
spectra, muons with large transverse momenta, cosmic ray 
composition and cosmic ray anisotropics on multiple an- 
gular scales. The cosmic ray anisotropics, the first time 
measured in the Southern sky, are drawing a lot of interest 
but have not yet found an explanation. 

IceCube can be used as a unique instrument to measure 
transient events, such as supemovae, GRBs and sun flares. 
This already led to results on heliosperic physics. 

Looking into the future: it seems as if the discovery of cos- 
mic high energy neutrinos might need some more years, in 
which the existing detectors will be exploited, improved 
and extended. The first, already accomplished, exten- 
sion was DeepCore opening the way to low energy phe- 
nomena such as neutrino oscillations, low mass WIMPs 
and supernova physics. A new low energy extension 
with very dense optical sensor instrumentation to allow for 
Cherenkov imaging in a megaton scale detector is studied, 
an interesting physics application being the search for pro- 
ton decay [51]. At the high energy end: radio and acoustic 
extensions are studied to reach the sensitivity for GZK neu- 
trino fluxes [52, 53] and to extend the air shower detection 
capabilities [54]. 

Acknowledgements : I would like to thank the many peo- 
ple in the IceCube Collaboration who helped me preparing 
the talk and the proceedings. 

References 

[1] The IceCube Collaboration, The IceCube Neutrino Ob- 
servatory I- VI: IceCube ICRC 2011 contributions: 
I: Point Source Searches, arXiv:l 11 1.2741; II: All 
Sky Searches: Atmospheric, Diffuse and EHE, 
arXiv:l 11 1.2736; III: Cosmic Rays, arXiv:l 11 1.2735; 
IV: Searches for Dark Matter and Exotic Par- 
ticles, arXiv:l 11 1.2738; V: Future Developments, 
arXiv:l 11 1.2742; VI: Neutrino Oscillations, Supernova 
Searches, Ice Properties, arXiv: 1111 .273 1 . 



[2] Felix Aharonian, "Probing cosmic ray accelerators 
with gamma rays and neutrinos", invited talk, these pro- 
ceedings. 

[3] IceCube Collab., ICRC contr. 807, arXiv:l 11 1.2735. 
[4] R. Abbasi et al. (IceCube), NIM A601 (2009) 294. 
[5] IceCube Collab., ICRC contr. 333, arXiv: 1111 .273 1 . 
[6] IceCube Collab., ICRC contr. 899, arXiv:l 11 1.2735. 
[7] IceCube Collab., ICRC contr. 1235, arXiv:ll 11.2741. 
[8] R. Abbasi et al. (IceCube), PRD 84 (2011) 082001. 
[9] IceCube Collab., ICRC contr. 329, arXiv: 1 1 1 1 .273 1 . 
[10] IceCube Collab., ICRC contr. 324, arXiv: 1111 .2736. 
[1 1] F. Kislat, Astrophys. Space Sci. Trans. 7 (201 1) 175. 
[12] R. Abbasi, et al. (IceCube), ApJ 732 (2011) 18. 
[13] IceCube Collab., ICRC contr. 784, arXiv:ll 11.2741. 
[14] R. Abbasi, et al. (IceCube), arXiv: 1104.0075. 
[15] E. Waxman. NP B Proc. SuppL, 118 (2003) 353. 
[16] D. Guetta et al., ApP 20 (2004) 429. 
[17] IceCube Collab., ICRC contr. 764, arXiv: 1111 .2741 . 
[18] R. Abbasi et al. (IceCube), PRE 106 (2011) 141101. 
[19] IceCube Collab., ICRC contr. 445, arXiv:ll 11.2741. 
[20] IceCube Collab., ICRC contr. 535, arXiv: 1111 .2741 . 
[21] IceCube Collab., ICRC contr. 334, arXiv:ll 11.2741. 
[22] IceCube Collab., ICRC contr. 736, arXiv: 1111 .2736. 
[23] Ahlers et al., ApP 34 (2010) 106. 
[24] E. Waxman and J. Bahcall, PRD 59 (1998) 023002. 
[25] IceCube Collab., ICRC contr. 759, arXiv: 1111 .2736. 
[26] IceCube Collab., ICRC contr. 1097, arXiv: 1111 .2736. 
[27] IceCube Collab., ICRC contr. 949, arXiv: 1111 .2736. 
[28] R. Abbasi, et al. (IceCube), PRD 83 (2011) 092003. 
[29] IceCube Collab., ICRC contr. 773, arXiv: 1111 .2736. 
[30] IceCube Collab., ICRC contr. 778, arXiv: 1111 .2736. 
[31] IceCube Collab., ICRC contr. 327, arXiv: 1111 .2738. 
[32] IceCube Collab., ICRC contr. 1 178, arXiv: 1111 .2738. 
[33] IceCube Collab., ICRC contr. 1024, arXiv:l 11 1.2738. 
[34] IceCube Collab., ICRC contr. 734, arXiv: 1111 .2738. 
[35] IceCube Collab., ICRC contr. 85, arXiv:l 11 1.2735. 
[36] J. R. Horandel, ApP 19 (2003) 193. 
[37] IceCube Collab., ICRC contr. 323, arXiv: 1111 .2735. 
[38] IceCube Collab., ICRC contr. 662, arXiv: 1111 .2735. 
[39] G. D. Barr et al., Phys. Rev. D70 (2004). 
[40] R. Abbasi et al. (IceCube), ApJ 718 (2010) L194. 
[41] R. Abbasi et al. (IceCube), ApJ 740 (2011)16. 
[42] IceCube Collab., ICRC contr. 305, arXiv:l 11 1.2735; 

Abbasi et al. (IceCube) arXiv: 1109. 1017, acc. by ApJ. 
[43] IceCube Collab., ICRC contr. 306, arXiv: 1111 .2735. 
[44] IceCube Collab.,ICRC contr. 308, arXiv:l 11 1.2735. 
[45] IceCube Collab., ICRC contr. 923, arXiv: 1111 .2735. 
[46] IceCube Collab., ICRC contr. 838, arXiv: 1111 .2735. 
[47] IceCube Collab., ICRC contr. 939, arXiv: 1111 .2735. 
[48] IceCube Collab., ICRC contr. 1 137, arXiv: 1111 .273 1 . 
[49] R. Abbasi et al. (IceCube), ApJ Lett. 689 (2008) 65. 
[50] IceCube Collab., ICRC contr. 921, arXiv:ll 11.2735. 
[51] IceCube Collab., ICRC contr. 325, arXiv: 1111 .2742. 
[52] IceCube Collab., ICRC contr. 1236, arXiv: 1111 .2742. 
[53] IceCube Collab., ICRC contr. 316, arXiv: 1111 .2742. 
[54] IceCube Collab., ICRC contr. 1 102, arXiv: 1111 .2742. 



12 



